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ABSTRACT 

We compare optical and high quality x-ray data for three bright ellipticals in the Virgo 
cluster, NGC 4472, 4649, and 4636. The distribution of total mass in NGC 4472 and 4649 
determined from x-ray data is sensitive to the stellar mass over a considerable range in galactic 
radius extending to r « r e , the effective radius. The agreement of x-ray and optically-determined 
stellar masses provides a unique verification of the stellar mass to light ratio which is essentially 
constant over the range 0.1 < r/r e < 1. However, for NGC 4636 the dark matter is important 
at all radii > 0.35r e . Evidently the dark to stellar mass ratio varies in quite different ways in 
ellipticals of comparable optical luminosity, implying that the radial structure of dark halos 
may not be universal. There is some evidence in NGC 4636 for additional support of the hot 
interstellar gas at r ^ 0.35r e ; either a field B ~ 10~ 4 G or a small (mechanically unstable) 
central region of high gas temperature (T ~ 10 K) is required. The global temperature 
structure in the hot interstellar medium of many recently observed ellipticals is very similar, 
reaching a maximum near 3 — 4r e . This feature, which may suggest a new structural scale in 
these galaxies, is inconsistent with current theoretical gas dynamical models. 



Subject headings: galaxies: elliptical and lenticular - galaxies: cooling flows - galaxies: 
individual(NGC4472, NGC4636, NGC4649) - xrays: galaxies 



INTRODUCTION 



Elliptical galaxies are large, bright stellar systems with considerable structural regularity. The stellar 
light profiles are well fit by the de Vaucouleurs law (Burkert 1993) and the stellar velocity dispersion, 
surface brightness, and effective radius r e are constrained to a fundamental plane (Djorgovski & Davis 
1987; Dressier et al. 1987). By contrast, the x-ray emission from hot interstellar gas in ellipticals exhibits 
enormous variations in L x /Lb (Eskridge et al 1995) and, to a lesser extent, in gas temperature (Davis & 
White 1996). 

We have re-examined optical and x-ray observations of ellipticals in preparation for a new series of gas 
dynamical studies of the evolution of the interstellar gas. Although current x-ray observations are sparse, 
we have found that they reveal more useful information about elliptical galaxies than is currently realized. 
Gas temperature profiles T(R) for elliptical galaxies recently determined with ROSAT and ASCA reveal a 
surprising uniformity for radii r < 10r e . The total mass and density of all gravitating matter, M tot (r) and 
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Ptot (r) , can be found from the variation of temperature and density in the hot interstellar gas by assuming 
that the gas is in hydrostatic equilibrium in the galactic potential. 

We show here that the total mass determined from the x-ray gas for two bright Virgo ellipticals (NGC 
4472 and 4649) is in excellent agreement with a de Vaucouleurs profile for 0.1r e r < r e using optically 
determined mass to light ratios. Remarkably, in another bright Virgo elliptical NGC 4636 dark matter 
contributes substantially to the total mass within r e and the total mass found from x-ray observations is 
lower than the known stellar mass for r ^ 0.35r e , implying that thermal pressure is not the only support for 
the hot gas. These results suggest that AXAF x-ray observations of elliptical galaxies will provide powerful 
new constraints on both the interstellar physics and the nature of stellar populations in these galaxies. 

2. SCALE OF THERMAL STRUCTURE IN ELLIPTICALS 

Figure 1 shows a plot of the radial variation of gas temperature with projected galactic radius R/r e for 
six early type galaxies. The remarkable feature that all galaxies share is a positive temperature gradient 
out to about ~ 3r e followed by a leveling off or gradual decrease toward larger radii. Data for Figure 1 are 
taken from the following sources: NGC 1399: Jones et al. (1997); NGC 5044: David et al. (1994); NGC 
4636: Trinchicri et al. (1994); NGC 507: Kim & Fabbiano (1995); NGC 4472: Irwin & Sarazin (1996); NGC 
4649: Trinchicri et al. (1997). The temperature profile observed with the ASCA satellite by Mushotsky et 
al. (1994) for NGC 4636 is in excellent agreement with that of Trinchicri et al. (1994) for 20 percent solar 
abundance. Additional properties of these galaxies are listed in Table 1. It is remarkable that T(R) is so 
consistent in spite of the wide range of luminosities and cosmic environments among these galaxies. NGC 
4636, 4649 and 4472 arc in the Virgo cluster while NGC 1399, 5044, and 507 are the brightest galaxies in 
small groups or clusters. NGC 507 has an exponential surface brightness profile that is consistent with a 
very massive, face-on SO galaxy (Magrelli, et al. 1992). 

The temperature maximum observed in ellipticals at ~ 3r e suggests a new previously unrecognized 
structural scale length. This cooling of gas within 3 — 4r e is not a natural result of galactic cooling flows 
as many have suggested. In theoretical models, the radial gas velocity in the hot interstellar gas is very 
subsonic since otherwise the gas replenishment time would be impossibly short and L x would be much too 
low. Straightforward spherical models of the hot interstellar gas (e.g. Fig. 3 of Brighenti & Mathews 1996) 
typically result in subsonic inflow with dT/dr < for r > 0.1r e . Although the gas loses thermal energy by 
emitting the x-rays observed, it is heated by compression as it descends deeper into the galactic potential, 
i. e. dT/dr < 0. In this sense the term "galactic cooling flow" is a misnomer. In "mass drop-out" models, 
in which gas is removed from the inner flow, the remaining gas is actually heated since the Pdv work 
done by gravity must raise its temperature even higher to provide hydrostatic support for gas at larger 
radii. However, this rise in T(r) is fully compesated when radiation from the cooler, dropping-out gas 
is included (e.g. White & Sarazin 1987). In any case, mass dropout is usually invoked at small galactic 
radii, less than the scale of dT(r)/dr < in Figure 1. The gas temperature is also sensitive to heating by 
supernova explosions, but this heating should be proportional to the stellar density so dT/dr < at r < r e 
would again be expected. To our knowledge no theoretical cooling flow model can produce temperature 
profiles resembling those in Figure 1. Fortunately, regardless of the possibly complex thermal history of the 
interstellar gas, the observed variation of T and p can be used to determine the underlying mass distribution 
as we now discuss. 
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3. COMPARISON OF X-RAY AND OPTICAL DATA FOR THREE VIRGO GALAXIES 

NGC 4472, 4649 and 4636 listed in Table 1 are all members of the Virgo cluster. We adopt a distance 
of D = 17 Mpc for all three galaxies, similar to the Malmquist-corrected value D = 17.2 ± 1.9 Mpc of 
Gonzalez & Faber (1997). For any particular Virgo galaxy the distance is uncertain by ~ 20 percent due to 
the unknown location of the galaxy along the line of sight through the cluster. 

The total mass M tot (r) in these galaxies can be found from the condition for hydrostatic equilibrium: 



where the last term allows for the possibility of magnetic pressure P m = B 2 /8ir and 4> m — P m / P is the ratio 
of magnetic to gas pressure. The proton mass is m p and we assume p = 0.63. The strong negative density 
gradient d\ogp/d\ogr is expected to be the largest of the three derivatives. Only the shape of p(r), not its 
absolute normalization, influences M tot (r). The density distribution p(r) in the hot interstellar gas can be 
found from the x-ray surface brightness distribution. T(r) and possibly also <f> m (r) must be known to some 
precision in order to determine reliable total masses. The total mass density is ptot(r) = dM tot /A r Kr 2 dr. 

Density profiles are available for all three ellipticals from Einstein HRI observations (Trinchicri, 
Fabbiano & Canizares 1986) and from ROSAT PSPC (Trinchicri et al. 1994; Irwin & Sarazin 1996; 
Trinchieri, Fabbiano, & Kim 1997). Fortunately there is a substantial range of angular scale over which these 
data sets overlap so values of p(r) from Einstein and ROSAT data can be renormalized to agree. Having done 
this we fit n(r) — p(r)/m p with a sum of functions n(r) = En^r) where n^(r) = n a (i)[l + (r/r (i)) p W] _1 
and the temperature is fit with T(r) — 2T m [r m /(r + r ot ) + (r/r m ) 9 ] -1 . The temperature T(R) observed at 
any projected radius R is an average along the line of sight weighted by p 2 and differs in principle from the 
temperature T(r) at physical radius r. However, p 2 is a very steep function of galactic radius and we find 
T(r) « T(R) within 10 percent, sufficient for our purposes here considering the observational uncertainties 
involved (Fig. 1). Figures 2a, 2d and 2g show our fits to the x-ray data. (The parameters for these fits are: 
NGC 4472: n a {i) = 0.095, 0.00597, -0.0004; r (i) = 0.107, 0.95, 10.; p(i) = 2.0, 1.14, 1.19; T m = 0.75; r m = 
0.5; r ot = 0.75; q = 0; NGC 4636: n Q = 0.151; r Q = 0.172; p = 1.57; T m = 0.5375; r m = 0.475; r ot = 0.6; 
q = 0.0; NGC 4649: n a (i) = .1, .0014; r a (i) = .15, 3.0; p(i) = 1.8, 3.0; T m = 0.9; r m = 4.0; r ot = 4.0; q = 
0.0 with radii in r e , densities in cm -3 and temperatures in 10 7 K.) Total masses M tot (r) (Figs. 2b, 2e and 
2h) and the corresponding total mass densities ptot{r) (Figs. 2c, 2f, and 2i) are determined with 4> m (r) = 0. 
Our values for M tot (r) in Figure 2b are in satisfactory agreement with those of Irwin and Sarazin (1996) 
(scaled to D = 17 Mpc) which are based only on ROSAT data; our M tot (r) in Figure 2h is in excellent 
agreement with Mushotsky et al. (1994). 

It is of particular interest to compare M tot (r) and ptotir) with corresponding stellar values. The total 
mass M* t is found from Lb using "stellar" mass to light ratios of van der Marel (1991), determined by 
comparing stellar velocities measured out to r a i, s « 0.5r e with solutions of the 2D Jeans equations in slowly 
rotating model galaxies. This "stellar" M* t /L B is sensitive to all mass within r obs , and may contain a 
component of dark matter which can be checked a posteriori (see below). In general mass to light values 
determined only in the core region (e.g. Faber et al. 1997) are larger; central black holes may account for 
this increase but typical holes may not be sufficiently massive (Kormendy, J. & Richstone, D. 1995; Faber 
et al. 1997). M*(r) and p*(r) in Figure 2 are evaluated with de Vaucouleurs profiles (Young 1976). 

The Structure of NGC 44^ an d 4649: We find it quite remarkable that the total mass and density 
indicated by the x-ray observations agree almost exactly with stellar values in the range 0.1 < r < lr e 
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where the observational data is of high quality (Figs 2b and 2e) . Discounting a conspiracy of compensating 
errors, this superb agreement is possible only (i) if the stellar mass to light values determined by van der 
Marel (1991) are essentially correct and constant for stars out to ~ lr e , and (ii) if magnetic pressure or 
rotation contributes little or nothing to the total pressure support of gas in 0.1 < r < lr e . The total masses 
of dark and stellar matter are equal at ~ 2.5r e which includes r oos so van der Marel's M*/Lb are likely to 
be totally stellar. Dark matter dominates at r/r e > 2.5. At r = 10r e the total mass to light ratio is ~ 78 
(NGC 4472) and - 110 (NGC 4649) and M tot /L B is undoubtedly higher at r/r e > 10. The stellar to dark 
halo transition is particularly striking in the density plots (Figs 2c and 2f). None of the qualitative features 
in Figure 2 are changed if (M*/ 'LB)core values are used from Faber et al. (1997), but the overall fit in 
0.1 < r/r e < 1 is of lower quality. The discrepancy M to t < M* at r < 0.1r e may result from observational 
inaccuracies in this region although M tot (r) would be underestimated if a large magnetic field were present. 

The Structure of NGC 4636: Our total mass M tot (r) for NGC 4636 (Fig. 2h) is in excellent agreement 
with the mass determined by Mushotsky et al (1994) who also assumed D = 17 Mpc. A de Vaucouleurs 
profile is used to determine M*(r) and p*(r); we have verified that NGC 4636 satisfies a de Vaucouleurs 
profile by plotting the surface brightness data of Peletier et al (1990) against i? 1 / 4 . 

The most surprising result evident from Figures 2h and 2i is the relative dominance of dark matter 
in NGC 4636 as compared to NGC 4472 and 4649. The dark mass M dark = M tot - M* in NGC 4636 
becomes equal to the stellar mass M* at about r = 1.2r e and Mdark ~ 0.83M* at r = r e where the total 
mass to light ratio is 19. This large amount of dark matter may be unusual since it has been difficult until 
recently to find evidence of dark matter in ellipticals from stellar velocities measured within r e (e.g. Carollo 
et al. 1995). The "stellar" mass to light ratio M^/Lb from van der Marel is determined from velocity 
observations only out to r oos /r e — 0.45, where M* ^> Mdark (Fig. 2h). It is remarkable that there is no 
change in slope in Figure 2h as M tot (r) crosses below M*(r) at r k 0.33r e since all observations should 
be reliable at this radius. This failure of the x-ray determined mass to detect the stellar mass suggests 
large magnetic fields, unusually high gas temperatures or rotational support in this region. Magnetic fields 
B~9x 10~ 5 G would be required at r = 0.1r e to support the hot gas against the stellar mass within (see 
Mathews & Brighcnti 1997). Or a curious central temperature inversion may exist: a mean temperature 
within 1 kpc (0.20') of 0.9 x 10 7 K is implied. This inversion is constrained so that the mean projected gas 
temperature within 1', 0.75 x 10 7 , is within the observed limits 0.667 + 0.029/ - 0.041 x 10 7 K. However, 
the high temperature interpretation for M to t < M* in NGC 4636 is unlikely since it is buoyantly unstable; 
dT/dr is supcradiabatic for r 0.08r e « 0.7 kpc. Alternatively, M*/Lb could decrease with galactic radius 
so that M* is lower than we think near r/r e w 0.1 in Figure 2h, but this seems unlikely in view of the 
uniform M*/L B implied by Figures 2b and 2e. None of these conclusions are changed if the core value 
(M*/ L B ) core = 12.69 is used instead. The different distribution of dark matter at r < r e in NGC 4636 may 
provide evidence against a universal dark halo structure as proposed by Navarro, Frenk & White (1996), 
but for r >r e the dark halos are all very similar having slope pdark °c r~ 19 at r <~ 10r e . The total mass to 
light ratio near the outer limit of the x-ray observations shown in Figure 2g, r = 12.5r e , is about 125. 

How do these Virgo ellipticals differ in other respects? NGC 4636 is a relatively isolated galaxy quite 
far (> 3 Mpc) from the core of Virgo. Nevertheless, its x-ray image shows some azimuthal asymmetry at 
r > 5r e , but this emission may arise from a different source (Trinchieri et al. 1994). It also has a broad, 
faint stellar distribution characteristic of CD galaxies. NGC 4472 is optically the brightest galaxy in the 
Virgo cluster, lies within a small subgroup, and is interacting with a nearby dwarf irregular galaxy, UGC 
7636. At radii > 3r e the x-ray isophotes show a tail-like structure which may result from the motion of 
NGC 4772 through the Virgo cluster medium (Forman et al 1985; Irwin and Sarazin 1996). NGC 4472 has 
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a small kinematically decoupled (non-rotating) core (Davies 1989) within <~ 0.096r e . NGC 4649 appears to 
be close to the center of Virgo. However, none of these galaxies shows evidence of unusual x-ray or optical 
structure at 0.1 < r/r e < 3 where the difference in dark to stellar mass is so apparent in Figure 2. 

4. CONCLUSIONS 

(1) A previously unrecognized spatial scale at ~ 3r e is suggested by the similarity of gas temperature 
profiles in many recently observed ellipticals. (2) The stellar mass to light ratios for NGC 4472 and 
4649 have been verified by x-ray observations; the interstellar gas temperature is also correct. (3) The 
mass to light ratio in NGC 4472 and 4649 determined from stellar motions are constant over the range 
0.1r e < r < r e . (4) In NGC 4472 and 4649 the hot interstellar gas is supported out to r ~ r e by thermal gas 
pressure, other means of support such as magnetic pressure or rotation are not evident. (5) Hot interstellar 
gas in the center of NGC 4636 may be supported by magnetic stresses (B ~ 10~ 4 G), a transient unstable 
region of unexpectedly high temperatures (T ~ 10 7 K), or rotation. (6) Elliptical galaxies with comparable 
luminosities can have quite different distributions of dark matter relative to stellar matter. (7) In NGC 
4636 the mass of dark matter is comparable to that in stars at r <~ r e ; in NGC 4472 and 4649 the mass of 
dark halo matter is negligible at r ~ r e . (8) The "stellar" mass to light ratios for these three ellipticals 
as determined by van der Marel (1991) are unlikely to be contaminated by dark matter from the galactic 
halos. AXAF observations will be useful in interpreting the variation of M^/Lb along the fundamental 
plane (Pahre & Djorgovski 1997) and in verifying central magnetically supported regions. 
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Fig. 1. — Radial dependence of gas temperature with projected radius in five early type galaxies all showing 
maxima at ~ 3r e . 

Fig. 2. — (a,d,g) Observations and our fits to gas temperature (top) and density (bottom) for Einstein HRI 
(filled circles) and ROSAT (open circles) data; (b,e,h) solid curves: M to t', points are values of M tot from 
Irwin and Sarazin (1996) for 4472 and Mushotsky et al. (1994) for 4636; long-dashed curves: stellar mass 
M„; short-dashed curves: dark mass Mdark = M to t — M* which loses accuracy when > M tot ; dot-dashed 
curves: total mass of hot gas; (c,f,i) solid curves: total mass density p to t', long-dashed curves: stellar mass 
density p* with a horizontal cut at the core or break radius r^; short-dashed curves: dark mass density. 
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Tablc 1. OPTICAL AND X-RAY PROPERTIES OF SIX ELLIPTICALS 
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a Type from RC3. 

b Distanccs NGC 507 & 5044 from Faber et al. (1989); NGC 1399 from Faber et al. (1997); 
NGC 4636 and 4472 from Gonzalez & Faber (1997). 

c Scaled to the distance in the second column. 
d Mass to light ratios from van der Marel (1991). 
e Break radii into core region from Faber et al.(1997). 
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